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X-ray powder photographs have been made of an electrically heated nickel ribbon sealed 
into a thin-walled glass bulb. X-ray measurements extending from 450 to 1200°C showed no 
new high temperature crystal form of nickel. The thermal expansion of the face-centered cubic 
lattice, present throughout the range investigated, is in agreement with macroscopic thermal 


expansion data. 


F the periodic group of elements, iron, 
nickel, cobalt, both iron' and cobalt? have 
two high temperature crystal forms in addition 
to the one usually occurring at ordinary tempera- 
tures. It would seem likely therefore that nickel 
would have another crystal form in addition to 
the usual face-centered cubic one. A hexagonal 
close-packed lattice has been reported by G. P. 
Thomson?’ as the result of electron diffraction 
measurements On a thin film of nickel sputtered 
on rocksalt in an atmosphere of argon. Attempts 
to obtain this form, or any other than the cubic 
form, by rapid quenching from high temperatures 
of nickel bars have not been successful. Such 
negative results in quenched bars do not elimi- 
nate the possibility of a high temperature crystal 
form since such a form may well revert to the low 
temperature form during even the most rapid 
quenching. Such is the case with pure iron. 
Recently a simple technique has been de- 
veloped by the writer whereby x-ray powder 
diffraction photographs of specimens in a vacuum 
or in an inert gas may readily be obtained. It 
seemed worth while therefore to make some 


1A. Westgren and A. E. Lindh, Zeits. f. physik. Chemie 
98, 181 (1921). 
 *S. B. Hendricks, M. E. Jefferson and J. F. Shultz, 
Zeits. f. Kristallographie 73, 376 (1930). 

°G, P. Thomson, Proc. Roy. Soc. A125, 358 (1929). 


measurements of a pure nickel ribbon at varying 
high temperatures to determine whether a second 
crystal form could be found and, in any event, to 
determine the dimensions of the crystal lattice at 
such temperatures. It seemed also that the 
method might be of value for other measurements 
of a similar character in case it proved itself 
capable of giving satisfactory results. 


EXPERIMENTAL PROCEDURE 


The x-ray measurements were made on a 
ribbon of very pure nickel bent into a hairpin 
filament and sealed into a thin-walled lime glass 
bulb filled with argon. The bulb was 2.5 centi- 
meters in diameter and the walls were less than 
0.01 cm thick. It was found that quite satis- 
factory powder photographs with copper radi- 
ation could be taken of a metal ribbon enclosed in 
such a bulb. The nickel used was formed by 
reducing pure nickel oxide in hydrogen, the 
resulting powder being then melted into an ingot 
in a vacuum furnace. This ingot was then 
swaged and drawn into a wire 10 mils (0.025 cm) 
in diameter and the latter was finally rolled out 
into a ribbon of dimensions 1344 mils (0.004 
0.11 cm). All annealing during the process of 
drawing was done in a vacuum furnace. A 
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spectroscopic analysis of the nickel showed only a 
trace of cobalt and a trace of silicon as impurities. 

The nickel filament was mounted in the thin- 
walled bulb and the latter was baked, after 
which the filament was degassed by heating and 
the bulb then filled with argon and sealed off. 
The x-ray measurements were made with the 
bulb mounted on an ordinary crystal spectrome- 
ter which was fitted with a precision camera of 7 
cm diameter. The film was placed in the camera 
so that the entering beam passed through it as in 
the method of van Arkel.4 The nickel filament 
was heated by passing a current through it froma 
heavy duty storage battery, the conditions of 
current and voltage being kept quite constant 
during each run. About a dozen photographs 
were taken of the glowing filament at tempera- 
tures ranging from 450 to 1200°C. The x-ray 
beam through the camera slits was directed at 
one arm of the hairpin filament and this arm was 
carefully centered in the camera by the use of a 
short focus telescope after the proper filament 
temperature had been attained for any given run. 
The bulb and filament were oscillated by means 
of the usual rocking cam, and care was taken that 
the second arm of the filament did not enter the 
path of the beam at any time during the oscil- 
lation to give spurious lines. The partial blocking 
off from the film of a small section of the diffrac- 
ted rays by the second arm of the filament was of 
small consequence. 

For filament temperatures above 800°C, tem- 
perature measurements were made with an 
optical pyrometer. Below 800°C, the tempera- 
tures were only approximately determined. 
Pyrometer temperatures were converted into 
true temperatures by making the proper cor- 
rections for the thermal emissivity of a bright 
nickel surface. The relation between brightness 
temperature (}=0.665u) and true temperature 
from the International Critical Tables was used.* 
No correction was made for the absorption due to 
the thin-walled bulb since experiment showed 
such a correction to be negligible. Although it is 
fully realized that pyrometer measurements of 
other than black bodies may be subject to large 
errors, it is believed that with the precautions 


4A. E. van Arkel, Zeits. f. Kristallographie 67, 235 
(1928). 
5 Int. Crit. Tab. 5, 245. 
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taken in this case the final temperatures obtained 
are very close to the true temperatures. 


RESULTS 


Measurements at room temperature 


Widely differing values for the size of the unit 
cube for nickel at room temperatures are quoted 
from the lfterature in the Strukturberichte® of 
Ewald and Hermann. Such values for the length 
of the cube edge vary from 3.499A to 3.536A. It 
seems therefore worth while to record the result 
of a series of precision measurements made at 
room temperature with the very pure specimen of 
nickel available. After small corrections‘ were 
made for faulty centering and any possible 
change in the length of the film during develop- 
ment, the values obtained from three different 
specimens for the length of the cube edge of 
nickel at 20°C were, respectively, 3.517 9A, 3.518,A, 
3.5172A. The mean 3.518+0.001A is in good 
agreement with the value 3.519+0.001A chosen 
in the Strukturberichte as the most probable value 
from the data available. 


High temperature measurements 


All the x-ray photographs of the hot nickel 
ribbon throughout the range of temperatures 
investigated (450 to 1200°C) gave lines corre- 
sponding to a face-centered cubic lattice. No 
lines corresponding to any other possible modi- 
fication were observed on any of the photographs. 
One is therefore forced to the conclusion that, 
unlike iron and cobalt, nickel shows no structural 
lattice transformation within the region investi- 
gated, and that a face-centered cubic structure is 
characteristic of the metal throughout this 
temperature range. It seems probable therefore 
that the hexagonal crystal form observed by G. 
P. Thomson’ cannot be obtained from thermal 
effects alone but probably is dependent for its 
existence upon the conditions of deposition of the 
sputtered film and possibly upon the base 
material upon which the nickel is sputtered. 

As is to be expected, due to thermal expansion, 
the dimensions of the face-centered cubic nickel 
lattice were found to increase with temperature. 
The values obtained for the length of the cube 
edge of the observed face-centered structure are 


6 Ewald and Hermann, Strukturberichte. 
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CUBE EDGE OF NICKEL LATTICE IN ANGSTROMS 


Fic. 1. Curve showing expansion of unit cube of nickel 
lattice with temperature. Full line from direct x-ray 
measurements, dotted line from calculations from thermal 
expansion data. 


plotted in Fig. 1 against the corresponding 
corrected pyrometer temperatures between the 
temperatures 850 and 1200°C. Data from two 
different specimens are indicated by A and ©, 
respectively. The resulting curve is indicated by 
the full line in the diagram. Average values 
obtained from the curve for the length of the 
cube edge at various temperatures are tabulated 
for convenience in Table I. 

A comparison of the values obtained above 
with values computed from the usual thermal 
expansion data seems desirable since, for single 
crystals of bismuth,’ discrepancies have been 
found between the macroscopic and the micro- 
scopic expansion data. No accurate results for 
the thermal expansion of nickel could be found 
covering the complete range 850 to 1200°C. 
Hidnert,* however, has carefully determined 


as 2 Goetz and R. C. Hergenrother, Phys. Rev. 40, 643 
). 
Peter Hidnert, Bur. Standards J. Research, Tech. 
Paper 257, 1306 (1930). 


TABLE I. Values of length of cube edge of nickel lattice at 
measured temperatures. 


Length of cube 
edge calculated 
from thermal 


Length of cube 


Temperatures edge of nickel 


degrees C lattice in A expansion data 
20 3.518+0.001 
900 3.568 +0.002 3.568 
1000 3.575 +0.002 
1100 3.583 +0.002 
1200 3.592 +0.002 


average expansion coefficients for a specimen of 
99.94 percent pure nickel for a number of 
temperature ranges between the limits 25 and 
900°C. Such data will serve for a calculation of 
the cube edge of the nickel lattice at tempera- 
tures up to 900°C, if we assume that the lattice 
expansion is exactly the same as the large scale 
expansion. The dotted curve in Fig. 1 is derived 
thus from Hidnert’s data where the value of the 
cube edge at 20°C is taken as 3.518A. This curve 
is seen to be a smooth continuation of the full line 
curve based on direct x-ray measurements. The 
value of the length of the cube edge at 900°C 
calculated from Hidnert’s average expansion 
coefficient for the full range 25—900°C is found to 
be 3.568A. (Table I.) This is just the value 
determined by x-ray measurements. It would 
seem therefore that within the limits of error of 
the measurements, the average expansion of the 
nickel lattice from 20 to 900°C as determined by 
direct x-ray measurements is the same as that 
calculated from macroscopic expansion data. 
The writer wishes to extend his thanks to Mr. 
E. F. Hennelly of the General Electric Research 
Laboratory for furnishing the samples of nickel 
used and to Mr. E. T. Asp for assistance in 
obtaining some of the x-ray photographs. 
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The Charges on Droplets Produced by the Spraying of Liquids as Revealed by the 
Millikan Oil Drop Method 


SEVILLE CHAPMAN,! University of California 
(Received February 24, 1934) 


HE electrification of liquids by spraying is a 

well-known phenomenon which has im- 
portant industrial as well as purely scientific 
interest. One phase of the problem, the investi- 
gation of spray electrification of water, has been 
thoroughly studied by P. Lenard? and his pupils. 
His researches indicated that the phenomena of 
cataphoresis of gas bubbles, waterfall electricity, 
electrification from jets and atomization by 
impact and bubbling are all intimately connected 
with the surface conditions of liquids, involving 
the so-called electrical double layer of the 
surface, and any electrolytic ions contained 
therein. 

According to Lenard and also to Coehn’ this 
stratification is determined by surface tension 
and dielectric constant of the liquid. Accordingly, 
while water should show a most pronounced 
effect, substances of lower dielectric constant 
should also show the effect, though to a lesser 
degree. Studies of cataphoresis of gas bubbles 
seemed to indicate that the effect, while present 
in water, was not observable in liquids of lower 
dielectric constant. Dyk* found, in such liquids, 
that before marked cataphoresis of air bubbles 
occurred, the turbulence caused by ionic move- 
ments in the high fields set in, thereby making 
further study impossible. Spray electrification, 
however, is daily observed, even in such sub- 
stances as pure hydrocarbons. Lenard and his 
pupils found that for water, the small droplets 
were negatively charged, while the larger drops 
were positively charged if ions were present; 
otherwise, they were neutral. In general, Biihl’ 
confirms these findings, but Busse® found that 
some positively charged drops were produced in 

1 Introduced by Leonard B. Loeb. 

2 Lenard, Ann, d. Physik 47, 463-524 (1915). 

§’Coehn, Ann. d. Physik 43, 1048-1078 (1914). 

‘ Dyk, Phys. Rev. 31, 913 (1928). 


5 Biihl, Koll. Zeits. 59, 346-353 (1932). 
® Busse, Ann, d. Physik 76, 493-533 (1925). 


pure water, although not in such quantity as the 
negative drops. 
_ As a quantitative method for a study of this 
spray electrification, it seemed that the Millikan 
oil drop method would give valuable information 
on the effect of drop size as related to charge. 
For the purpose of this study, a conventional 
Millikan oil drop apparatus was employed, the 
droplets entering the condenser space by means 
of gravity through a hole in the upper plate. 
The condenser plates were separated by 0.8 cm. 
A potential from 2000 volts downward, could be 
applied to the plates. Light from a 1000 watt 
projection lamp, filtered with a CuCl solution in 
a water cell, was used for illumination. The 
telescope had a 10 cm focal length and a mag- 
nification of 20 times. The droplets were pro- 
duced by a brass atomizer of a conventional 
design, in which the liquid was drawn up through 
a fine brass tube and atomized by a high velocity 
stream of clean, filtered, compressed air issuing 
from a nozzle at right angles, pressures of 6-8 
Ibs./in.? being used. The substances used, to- 
gether with the purification resorted to, are 
listed below: 


1. Ameroil, a pure, inert medicinal hydrocarbon, largely 
octane; taken from drugstore samples. 

2. Nitrobenzene; highly purified for use in Kerr cells, as 
prepared by H. J. White in his study of Kerr cells. This was 
especially pure nitrobenzene, which was exceptionally free 
from ions. 

3. Aniline; twice redistilled in air, using only the sharp 
middle fraction. 

4. Glycerine; the water heated out of it and distilled 
in vacuum, only the middle fraction being used. 

5. Ordinary distilled water. 


Further substances were not tried, as it was 
difficult to find any of significantly different 
properties which would not evaporate too 
rapidly. 

All measurements were made at atmospheric 
pressure and 20°C. With water, nitrobenzene and 
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Fic. 1. Charge as a function of drop size in nitrobenzene. 


aniline, considerable difficulty was encountered, 
due to the evaporation of the drops during 
measurements. This was only partially overcome 
by saturating the air inside the chamber with the 
vapor of the liquid. 

The results obtained in the case of nitrobenzene 
are summarized in the graph of Fig. 1. In this 
graph, the drops, characterized as positive or 
negative by a + or — sign, are plotted with 
their charge in electrons against the radius of the 
drop in centimeters. It must be noted that, 
owing to evaporation, accurate values of the 
charge in exact multiples of the charge of the 
electron could not be obtained. Consequently, 
the values were computed to the nearest whole 
multiple of e. The values are good to about 15 
percent. 

It is seen at once that the charge varies in 
magnitude with drop size, being greater for the 
larger drops. This holds for all substances, except 
Ameroil, where all charges are uniformly low. 
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It is to be emphasized that in this range of sizes 
no substance (including Ameroil) showed any 
correlation of drop size with sign of charge. 
Different substances, however, varied widely in 
the magnitude of the charge for approximately 
the same size of drops. The relation of the 
magnitude of the charge to dielectric constant 
and other properties of the substances are shown 
in Table I. For each substance in the table are 
given the dielectric constant, the surface tension, 
and the viscosity taken from standard tables; the 
radius of the smallest drops with their average 
charge, and the radius of the largest drops with 
their average charge. 

The drop sizes were limited to relatively 
narrow ranges of size. This is due to a number of 
causes. The mass of the drops increases as the 
cube of their radii. Thus, for larger drops the 
short time of fall under gravity, and the necessity 
of sufficiently high charges to make them alter 
their speed, precluded measurements of very 
large drops. This was especially important in the 
case of Ameroil, where the charges per drop are 


small, but was of little significance in the cases of 


the polar liquids, which gave more highly charged 
drops. Measurements on small drops become 
impossible, owing to the violent Brownian, and 
slight convectional movements, and in the case of 
the more volatile liquids to the rapid evaporation 
of the drops. This latter factor made the results in 
water yield orders of magnitude only. These 
measurements were, however, important in 
showing that in this range of drop size, water 
behaved in a fashion similar to the other 
substances. It is very possible that the actual 
magnitudes of the charges on the water drops 
could have been greater than those on nitro- 
benzene, although this is not certain. 


TABLE I. Charge on drops as related to other liquid properties. 


Smallest Drops Largest Drops 


Surface Viscosity Average Average 

Dielectric tension in in dyne Radius charge in Radius charge in 

Liquid constant dynes/cm sec. /cm? in cm electrons incm electrons 

1074 

Ameroil 1.9 27.0 0.054 4.01075 2 1 2 
Nitrobenzene 35.7 46.0 0.021 1.010~4 150 4 800 
Aniline 7.3 44.0 0.044 1.5x10-* 100 5 450 
Glycerine 56.2 65.2 0.083 5.0107 40 25 500 
Water 80.0 72.8 0.010 1.5x10-* 125 5 600 
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From the table it is seen that in the range of 
drop size used there is: 


1. No preference for either sign of charge. 

2. The larger drops of both signs of charge have the 
larger charges, and that the charge varies roughly in a 
linear fashion with drop size, except as previously noted in 
the case of Ameroil. 

3. That there is a rough parallelism between dielectric 
constant and charge magnitude for drops of approximately 
the same size. 

4, That although there appears to be no correlation with 
viscosity, a parallelism between drop size and surface 
tension exists which is, however, difficult to differentiate 
from that with dielectric constant. 


There was only one substance, Ameroil, which 
showed any uncharged droplets, and these were 
of all sizes. They formed about 25 percent of the 
drops observed. 

It is quite clear that the phenomena involved 
here differ radically from those observed by 
Lenard and his group. It may be that the 
mechanism of charged drop production is differ- 
ent for drops of the two ranges of size. Or, 
perhaps, in the present work the charges were 
due to the presence of ionized particles in the 
organic substances used. Lenard’s group sprayed 
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water and estimated drop size from the mobilities 
of the ions found in the air carrying the spray 
particles. They could not possibly have observed 
the slow droplets studied here. What they 
studied were either initially created minute drops 
(of about 10~? cm radius), or else larger droplets 
that had evaporated to the minute dimensions at 
which the charge on the droplet caused a balance 
between evaporation and condensation in the air 
blast used in the measurement. In the work on 
small drops a marked difference in mobility 
between positively and negatively charged drops 
was observed, and positively charged drops were 
found only in the case of salts dissolved in the 
water.’ 

It is clear that the relation between the 
charged droplets studied by Lenard and his 
group and those studied in this investigation can 
be determined only by a study of the mobilities 
and sign of the charges of drops of all sizes. To 
this end an ion mobility method of high resolving 
power, suitable to such an investigation, is being 
developed in this laboratory. 


7 On the contrary, Busse, reference 6. 
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Investigations of Turbulent Flow by Means of the Hot-Wire Anemometer 


F. L. WATTENDORF AND A. M. Kuetue, California Institute of Technology 
(Received March 24, 1934) 


I. INTRODUCTION 


URBULENT flow is characterized by an 

irregular fluctuating motion superimposed 
on a mean flow. Reynolds,! in his analysis of 
turbulent flow, expressed the shearing stresses 
for a parallel mean flow in the x-direction 


Try = —pu'd ; Tr: = —pu'w’, (1) 
where x’, v’ and w’ are the fluctuations in the 
velocity in the direction and perpendicular to 
the direction of the mean flow, respectively, 
p is the fluid density and the bar indicates that 
the mean value of the product over the time has 
been taken. It is seen, then, that the turbulent 
flow of a fluid at a point is characterized by the 
magnitude of the components of velocity fluctu- 
ation together with the correlation between 
them, and experimental evidence concerning 
these quantities is of great importance for the 
development of a theory of turbulent flow. 

The present. work gives measurements of a 
root-mean-square value of the fluctuation in the 
direction of the main flow, i.e., (u”)"/? in various 
types of turbulent flow. The method used was 
that of the compensated hot wire anemometer, 
the investigations were made for the following 
cases : 


1. Fully developed parallel flow. 

2. Boundary layer along a flat smooth wall. 

3. Fully developed curved flow. 

4. Working section of 10’ wind tunnel with and without 
artificial turbulence (in correlation with sphere drag tests 
and maximum lift measurements). 

5. Turbulence in the free atmosphere with airplane and 
blimp. 


The compensated hot-wire circuit used in the 
present work is based in principle on the circuit 


10. Reynolds, An Experimental Investigation of the 
Circumstances which Determine Whether the Motion of 
Water shall be Direct or Sinuous and of the Law of Resistance 
wa Channels, Phil, Trans. Roy. Soc. 174, 935 

). 


developed by Dryden and Kuethe** at the 
Bureau of Standards. The reader is referred to 
references 2 and 3 for details of the analysis, 
since only the results will be given here. 


II]. Discussion oF THEORY OF LAG or Hor 
WIRE AND COMPENSATION 


King’s* equation for the heat loss from an 
electrically heated wire in a steady stream of 
fluid is 

(2) 
where 1=current flowing through the wire, R 
=resistance of the wire, K =proportionality 
factor allowing for free convection, C=pro- 
portionality factor allowing for forced convec- 
tion, 7=temperature of the wire, 7) =tempera- 
ture of moving fluid. If now the heated wire is 
placed in a flow in which the speed fluctuates 
within small limits, it is found that the in- 
stantaneous temperature of the wire is a function 
both of the amplitude and the frequency of the 
fluctuations. Within narrow limits (about 15 
percent of the mean velocity) at a given fre- 
quency the amplitude of the temperature fluctu- 
ation is very nearly proportional to that of the 
speed fluctuation. If the amplitude remains 
constant and the frequency varies, as is shown 
in reference 2, the resistance fluctuation is cut 
down in the ratio 


1 jtan™ Moe 3 
(3) 


where .\/ is the time constant of the wire and is 
given by 7 
M=4.2ms(T—T (4) 


2H. L. Dryden and A. M. Kuethe, The Measurement of 
Wind Speed Fluctuations by Means of the Hot-Wire Ane- 
mometer, N. A. C, A. Technical Report 320 (1929). 

8H. L. Dryden and A. M. Kuethe, Effects of Turbulence 
in Wind Tunnel Measurements, N. A. C. A. Technical 
Report 342, (1930). 

*L. V. King, Phil. Trans, Roy. Soc. (London) A214, 
373 (1914), 
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Fic, 1. Schematic diagram of compensation circuit. 


In the above expressions w=27 times the fre- 
quency, m=the mass of the wire, 7=mean 
temperature of the wire, s=the specific heat of 
the material of the wire, Ro=resistance of the 
unheated wire. The phase lag expressed by the 
exponential term is important only when the 
actual shape of the fluctuations is being in- 
vestigated. 

The steady potential drop across the hot- 
wire terminals is balanced out by causing the 
wire to form one arm of a balanced Wheatstone 
bridge circuit. Then the fluctuations in voltage 
are impressed across the input of a vacuum tube 
amplifier. The root mean square of the current 
fluctuations in the output of the amplifier may 
be measured or the output may be used to 
actuate an oscillograph element. However, be- 
cause of the above-mentioned lag of the hot- 
wire, the results would not be a true measure of 
the velocity fluctuations except for very low 
frequencies. 

In order to counteract this lag, an electrical 
filter circuit of the type shown in Fig. 1 was 
included in the plate circuit of one of the amplifier 
stages. It consists of an inductance, a variable 
and a fixed resistance in series replacing the plate 
resistor. The drop in voltage across R, and L is 
passed on to the next stage. A complete dis- 
cussion of this circuit is given in reference 2. 
It will only be noted here that a measure of the 
theoretical effectiveness of the compensation is 
given by the term 


wL/(Rot+RitR2) (5) 
and that the attenuation factor is approximately 
R2/(Rot+Rit+R2). (6) 


It is seen that increasing the effectiveness of 


the compensation involves an increase in attenu- 
ation as well, and it becomes necessary then to 
adjust these two factors so that the performance 
of the entire arrangement will be as effective as 
possible without too much loss of amplification, 

In 1930 Ziegler® reported on an equivalent 
circuit which had been used with success at 
Delft. The compensation was obtained by the 
use of a condenser, a parallel and a series 
resistance. The arrangement was in the input 
circuit, the drop across the series resistance 
being passed on to the amplifier. Mock and 
Dryden® pointed out that if the compensation 
is to be placed between two of the amplifier 
stages the performance of the inductance circuit 
is superior to that of Ziegler’s arrangement. 

The fluctuations in resistance of a single hot- 
wire actually represent, to a first order, only u’, 
i.e., the velocity fluctuations in the direction of 
the main flow. This will be appreciated from the 
following analysis. Neglecting the lag of the wire, 
the instantaneous resistance corresponds to the 
velocity 


gw" 

where g is a factor, considerably less than unity, 
to allow for the fact that w’ fluctuations are 
parallel to the wire. Since the steady voltage 
drop across the wire is balanced out, we pass on 
to the amplifier an instantaneous voltage corre- 
sponding to the velocity a=q—g. The quantity 
measured by the mean square meter in the out- 
put corresponds to a? and 


since, by the rules for taking mean values 
gg=@. Now, expanding g by means of the 
binomial theorem, 


Squaring this quantity and subtracting from 


5M. Ziegler, The A pplication of the Hot-Wire Anemometer 
for the Investigation of the Turbulence of an Atr-stream, 
Proc. K. Akad. v. Wetenschappen Amsterdam 33.7 (1930). 

®W. C. Mock and H. L. Dryden, Improved Apparatus 
for the Measurement of Wind Speed Fluctuations, : oo 
A. Technical Report 448 (1932). 
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we get 
a=u", approximately. (7) 


Obviously, if the velocity fluctuations are a 
considerable percentage of the mean velocity at 
the point, more terms must be taken in the above 
expansion and we then measure an appreciable 
part of the v’ fluctuations. 


III. APPARATUS 


The hot-wire used at the Guggenheim Aero- 
nautics Laboratory, California Institute of Tech- 
nology (abbreviated GALCIT) was of platinum 
0.0147 millimeter in diameter and about three 
millimeters long. The diameter was measured by 
means of a travelling microscope and is probably 
correct to +0.0003 millimeter. The wire is 
welded across the terminals in the manner de- 
scribed in reference 3. The mounting is 3/32” 
drawn copper tube, bent at right angles one 
inch from the end nearest the hot-wire. An 
insulated wire running through the tube forms 
one terminal while a prong is fastened to the 
copper tube itself for the other terminal. Fig. 2 
shows one of the hot-wires mounted in a microm- 
eter screw for making traverses across a channel 
or through a boundary layer. 

Fig. 3 is a wiring diagram of the circuit used. 
The amplifier and part of the input circuit were 
built up at this laboratory, and the whole was 
adapted by the authors to the measurement of 
turbulent fluctuations. 


Fic. 2. Hot wire mounted on a micrometer traversing 
mechanism, 


The hot-wire in series with a one-ohm re- 
sistance and a milliammeter forms one arm of a 
Wheatstone bridge. The bridge is adjusted so 
that only the fluctuations in voltage across it 
are passed on to the input transformer of the 
amplifier. The wire is operated at a constant 
current measured by means of the potentiometer 
arrangement across the one-ohm resistance in 
series with the hot-wire. By means of this 
arrangement the current can be adjusted to 
+0.4 milliampere. 

The fluctuations are first amplified through 
two resistance-capacity coupled stages and then 
are fed into a three stage push-pull, transformer- 
coupled amplifier constructed by the General 
Radio Corporation. This unit is a.c. operated, 
the three stages using UY227, UX226 and 
UX250 tubes, respectively. The alternating 
current introduced a small hum in the amplifier 
but this was, at all times, small enough so that 
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Fic. 3. Wiring diagram of the GALCIT hot-wire apparatus. Since the'a.c. amplifier is of con- 
ventional type, the wiring diagram is not included. SC is a standard cell by means of which 
the voltage of E is determined. R,=10,000 ohms, L=2.45 henrys, R2=400 ohms (about), 
C=2 microfarads, A; is a Weston thermogalvanometer reading current squared. 
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its subtraction did not sensibly impair the 
accuracy of the results. With the compensation 
circuit removed the five stages gave an amplifica- 
tion which was constant with frequency to 
within five percent from 50 to 8000 cycles. 
This was considered a satisfactory framework 
in which to install compensation. 

The compensation circuit was designed with 
recourse to expressions (5) and (6), respectively. 
The constants of the circuit used were L=2.45 
henries, R»o=10,000 ohms, R;=10,000 ohms, 
R.=400 ohms (approximately). With these 
values (5) gives an error factor of 1.20«10-4 
from which we deduce that the error in com- 
pensation is one percent at a frequency of 185 
c.p.s., and seven percent at 500 c.p.s. The 
attenuation factor, calculated from (6) is 0.0196. 
For the scale at which measurements were made 
oscillograms showed that frequencies above 500 
c.p.s. were of secondary magnitude. 

Fig. 4 shows curves for the amplifier response 
with compensation, the lag of the hot-wire for 
M = 0.0036 and the result of the multiplication of 
these curves. The latter represents the overall 
response of the hot-wire, compensation and 
amplifier. Actually, the error factor mentioned 
above was only five percent at 500 c.p.s. and 
about ten percent at 800 c.p.s. 

The calibration procedure, giving the amplifier 
curve in Fig. 4, was as follows: The hot-wire was 
replaced by an equivalent constant resistance 
through which an alternating current of known 
frequency and amplitude, generated by a vacuum 
tube oscillator, was passed. The bridge was then 
connected across the input of the amplifier and 
the output was measured on a mean square 
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Fic. 4. Curves of hot-wire lag, gain of compensated am- 
plifier, and gain of hot-wire plus compensated amplifier. — 
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galvanometer. This was carried out for various 
frequencies from 50 to 1000 c.p.s., and for 
various values of the compensation resistance 
R: with and without inductance. By placing the 
same resistance across the terminals marked 
input and closing switch S,; the procedure was 
repeated. Then the ratio of these results to 
those obtained with the resistance in the bridge 
circuit (after applying a correction term re- 
sulting from the fact that the resistance was 
actually in parallel with the bridge circuit) gives 
the bridge factor, i.e., the factor by which the 
bridge cuts down the response because of the 
drop through the hot-wire, the one-ohm re- 
sistances and the milliammeter. It was found 
that the bridge factor did not vary by more than 
five percent in the entire range of frequencies, 
and from this it was concluded that the im- 
pedance of the input transformer was so large 
that it could be disregarded insofar as it might 
alter results. The value of the inductance used 
is really an effective inductance calculated from 
the slope of the amplification curve of the 
compensated amplifier and the amplification for 
L=0 and R.=1 ohm. The attenuator designated 
by R: has fifteen taps by means of which the 
amplification can be varied from full to prac- 
tically zero amplification. This was also cali- 
brated by means of the oscillator. 

The procedure of making turbulence measure- 
ments will be described briefly. The positions of 


the various switches can be determined from 


Fig. 3. The measurements were generally made 
at various positions along a lateral traverse of a 
channel, or through a boundary layer, and con- 
sisted of about ten measurements at as many 
stations. Before and after each traverse, it was 
necessary to measure the resistance of the wire 
at room temperature and to take a calibration of 
resistance of the wire against mean wind speed 
at the heating current used for the measure- 
ments. For the former, a very small current was 
passed through the bridge and by adjusting Ra, 
a balance was obtained. For the latter, the hot- 
wire was placed at a positien where the mean 
speed was known as a function of the static 
pressure at the wall of the channel. Maintaining 
a constant heating current, readings of the 
resistance of the hot-wire were taken at a series 
of wind speeds covering the range encountered 
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in the traverse. Then, plotting ?RRoa/(R— Ro) 
versus u'/? results in a straight line. That this 
is true is easily seen from a consideration of 
Eq. (8) below, which is Eq. (2) in a slightly 
different form. As is shown below, the slope of 
this line is used in the calculation of the velocity 
fluctuations. There was generally a slight diver- 
gence between the two calibration curves but the 
slopes were usually within 5 percent of each 
other. 

At each position of the traverse readings were 
taken of the mean resistance of the wire, the 
attenuator tap used and the root-mean-square 
value of the output as read by a thermal 
ammeter. Also at each position the time constant 
of the wire was calculated according to (4), and 
the compensation resistance adjusted to its cor- 
rect value so that the amplitudes of the fluctu- 
ations would be reproduced in the correct re- 
lation with each other. 

The reading of the thermal galvanometer 
cannot be interpreted directly in terms of ve- 
locity fluctuations since any voltage fluctuation 
across the wire is accompanied by a change in 
the heating current. Fig. 5 is a schematic 
diagram of the circuit considered. Since T— 7 
=(R—Ro)/Roaw Eq. (2) can be expressed in the 
form 


?RRoa/(R—Ro) =K+Cu}. (8) 


Differentiating this expression with respect to 7, 
R and u gives 


2iRRoa 
di— d 
R-Ry (R-—R,)? 


Cdu 


2u} 


and this, along with the relations 

E,(di/i?) = —(ri/rs+1)dR, 
dE=idR+Rdi, 
where 72=74+/;, gives the result 

G=2iRRya/(R—Ro), 


where 
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Fic. 5. The bridge circuit. 


In the above formula du is really [u’”|! since dE 
=(AE*)!. dE is found by dividing the output 
voltage by the overall amplification factor which 
includes the bridge factor previously discussed. 
The fact that the heating circuit was discon- 
nected for the calibration of the amplifier with 
the vacuum tube oscillator actually causes a 
very slight error. 


IV. DISTRIBUTION OF VELOCITY FLUCTUATIONS 
THROUGH A CHANNEL 


One of the fundamental problems in the study 
of turbulent flow is that of the mechanism of 
the transport of momentum in parallel or nearly 
parallel flow. Several theories have been ad- 
vanced for the solution of this problem, but for 
experimental checking it would be necessary to 
have measurements of the actual components of 
velocity fluctuations, u’ and v’. The present 
measurements give only the root-mean-square 
values of u’, the component parallel to the mean 
flow. Measurement of the perpendicular com- 
ponent v’ and the correlation between uw’ and v’ 
are in progress at the California Institute of 
Technology. 

Most of these measurements were made in a 
channel previously used by F. L. Wattendorf? 
for the investigation of curved flow. The channel, 
shown in Fig. 6, has a breadth of 5 centimeters 
and a depth of 90 centimeters so that the ratio 
of depth to breadth is 18:1. The ratio was 
chosen large in order to avoid, as much as 
possible, the disturbing influence of the top and 
bottom, or, in other words, so that the mean 
flow in the section at one-half the channel 
depth would be approximately two-dimensional. 


7 To be published at an early date. 
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CURVED CHANNEL 


Fic. 6. Diagram of curved channel showing the various 
measurement stations. 


The channel consisted of a bell-shaped en- 
trance, with honeycomb, leading into a straight 
section with parallel walls 285 centimeters long. 
The purpose of the straight section was to build 
up the flow into a fully developed turbulent 
state. This length is 57 times the channel 
breadth b, and previous investigations showed 
that the flow was in a fully developed state at 
the measuring section. The straight section, in 
turn, leads into a curved section, with the walls 
formed into concentric circular arcs, of inner 
radius 20 centimeters and outer radius 25 centi- 
meters. The curvature extends through about 
300° of arc, in order to obtain a fully developed 
curved flow. Fastened to the end of the curved 
section is a transition section leading to the 
propeller and motor unit. 

The measurement stations are also shown on 
the diagram. The hot-wire was projected through 
a hole in the tunnel wall, and its micrometer 
mounting (see Fig. 2) by means of which 
traverses were made, was securely fastened to 
the outside. A static orifice at each station 
measured the wall pressure at that point. 

Two different experimental techniques were 
used at section 9 which represents fully de- 
veloped straight flow. The values of u”/tn0x 
obtained from both methods are given in Fig. 7 
and it is seen that they check rather well with 
each other. The first method consisted of 
operating the hot-wire under conditions such 
that a constant value of M/, the time constant 
of the wire, would be preserved at all times. 
This meant that at every position during a 
traverse, the current through the wire had to be 
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Fic. 7. /t max and /tmax versus distance from the 
center for fully developed straight flow. 


adjusted to a value which would keep the 
ratio 7?/(R—Ro,) a constant. The second method 
involved the use of a constant heating current 
(usually 0.16 ampere) and calculating M/ and 
R2 for each point. The latter method was adopted 
as the standard procedure because of its greater 
ease of application. 


The distribution of u”/u?,.,.x for the fully 
developed straight flow is given in Fig. 7, and 
the significant thing is that a straight line can 
be passed approximately through all points 
except those near the wall and near the channel 
center. It is well known that for a fully developed 
parallel flow the friction varies linearly from a 
value of (b/2)(dp/dx) at the walls, where dp/dx 
is the pressure drop in the channel and 3 is 
the channel breadth, to zero at the center. 
The distribution of which, according to 
(1), is equal to —pu’v’/u*nax is shown also in 
Fig. 7. Comparison of the two linear distributions 
shows that 


approximately holds in the linear region.*® 


8 Since the completipn of this work an investigation by 
H. Reichardt (Zeits. {. Angew. Math. u. Mech. (13), 3, 
177-180 (1933)) has been reported, in which measurements 


of (u’)}/u are made in rectilinear fully developed flow at 
a low scale. Measurements at a low mean velocity with an 


uncompensated hot-wire are given. Plotting of —w’v'/w’ 
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The above results indicate a constant corre- 
lation between w’v’ and uw”. In the theory of 
turbulence the quantity —w’v’ /(w”)4(v”)! is used 
as a ‘‘correlation factor” and we need to know 
this quantity, along with the magnitude of the 
fluctuations, in order to fully characterize tur- 
bulent flow. If now we assume similarity of the 
turbulence mechanism, as proposed by von 
KArman,? which implies that the ratio of (w’)!/ 


Vey! / 


(v’”)} is constant, then u’v’/u’ is proportional to 
uv’ /(u”) Vy!” )}, and our results indicate therefore 
a constant correlation factor for straight flow. 
It may be mentioned that Prandtl’s’® mixing 
length theory contains the assumption that | u’| 
and |v’| are both proportional to /(du/dy), where 
| represents the so-called mixing length, so that 
this theory also implies constant correlation. 
A final proof of the assumption of constant 
correlation can only be found by direct measure- 


ment of (v”)! for which measurements are now 
being prepared in this laboratory. It may be 
mentioned that measurements by Fage," carried 
out by means of an ingenious optical method, 
show that for turbulent flow in a square channel, 
the maxima of wu’, v’, w’ are nearly equal in the 
region excluding the wall neighborhood. Also 
observations by Taylor’ of atmospheric tur- 
bulence, which represents a turbulent flow of 
large magnitude and low frequency, show that 
except near the ground w’ and v’ are nearly equal. 

For the purpose of investigating scale effect on 
the magnitude of the fluctuations, measurements 
were made at the channel center at section 9 for 
various mean speeds from 7 to 25 meters per 
second. Previous measurements of the mean 
velocity profiles showed the effect of Reynolds 
number on the velocity distribution to be small. 


The measurements of (”)}/umax, however, 
showed a definite decrease with increasing ve- 
locity, and the question is, on what physical 


over the section results in an approximately constant value 
of 0.24 over a large part of the region excluding the neigh- 
borhood of the walls and channel center. 

®Th. von Karman, Mechanical Similitude and Turbu- 
lence, Nachr. Ges. Wiss., Goettingen, 58-76 (1930). 

“L. Prandtl, Bericht tiber Untersuchungen zur Aus- 
gebildeten Turbulenz, Zeits. {. Angew. Math. u. Mech. 5, 
136 (1925). 

"A. Fage, Study of Turbulent Flow by Means of the 
Ultra- Microscope. Proc. Roy. Soc. A135, 828, 656 (1932), 

%G. I. Taylor, Turbulence, Roy. Met. Soc. J. 53, 201- 
212 (1927), 
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parameters is the scale effect based. Since the 
shearing stress is proportional to the product of 
the fluctuation components it would seem reason- 


able to assume a dependence of u” on the 
shearing stress at the wall ro. Accordingly, the 


curve of (u”)/tmax against R.N. was compared 
with that of 
Umax versus R.N. where K was arbitrarily chosen 
to bring the curves at the same level for purposes 
of comparison. Fig. 8 shows the result, and it is 
seen that the curves are similar, although the 
range of scale covered was not large enough to 
make any definite conclusions. 

However, it seemed possible that the magni- 
tude of the fluctuations in the channel center 
might also be dependent on entrance conditions. 
Accordingly, screens of various sizes were intro- 
duced just ahead of section 1, and the traverse 
at section 9 repeated. It was found that the 


distribution of (uv)! at section 9 was unaffected 
by the presence of the screens. It is well known 
that the mean velocity distribution in fully de- 
veloped turbulent flow is not affected by dis- 
turbances at the entrance of the channel. It 
follows, then, that if the fluctuation level were 
raised by the above procedure, it must be 
accompanied by a weakened correlation in order 
that the shearing stress remain unchanged. The 
results indicate, then, that in fully developed 
turbulent flow, the fluctuations in velocity have 
reached a statistical equilibrium and the corre- 
lation factor is unaffected by fluctuations intro- 
duced at the entrance. 

Traverses were also made at stations 2 and 3 
in the entrance section in order to study the 
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Fic. 8. The effect of scale on (u”)!/umax at the center 
of section 9, 
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Fic. 9. (u”))/umax and u“/tmax for two stations near 
the entrance of the channel. For comparison the results 
for fully developed straight flow are reproduced. 


development from boundary layer flow to fully 
developed turbulence. The results are shown in 
Fig. 9, along with the mean velocity profiles. 
Section 2 indicates a distinct boundary layer at 
each wall separated by a region of potential 
flow in the center. The magnitude of the fluctu- 
ations decreases rapidly through the boundary 
layer to a constant low value of about 0.5 percent 
in the region between the boundary layers. At 
section 3 the velocity profile still indicates two 
distinct boundary layers, separated by a rather 
narrow gap, but the rise in turbulence level at 
the center indicates that there is already con- 
siderable turbulent interchange between the 
boundary layers. It is significant that the magni- 
tude of the fluctuations near the wall is of the 
same order for all sections. 


\V. MEASUREMENTS IN THE CURVED CHANNEL 


Traverses of (wu)! as well as mean velocity 
were taken at various stations in the curved 
portion of the channel. The results are shown in 
Fig. 10. For convenience the ordinates are dis- 
placed vertically by the amounts indicated. 

The fluctuation profiles throughout the transi- 
tion region from fully developed straight to 
fully developed curved flow show two pre- 
dominating tendencies. The first is a general 
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flattening of the profile in the central region of 
the channel. This is closely analogous to the 
corresponding velocity profiles, which show that 
the mean flow in the corresponding region 
approaches the potential distribution ur=con- 
stant, where 7 is the radius at the considered 
point. The second tendency is that the gradient 


of (u)} is steeper at the outer wall and flatter 
at the inner wall than for straight flow. This is 
probably associated with the different states of 
stability at the two walls. At the outer wall the 
centripetal pressure gradient is directed away 
from the wall, giving rise to instability of the 
turbulent exchange, while at the inner wall 
the pressure gradient is directed towards the 
wall, which produces a stabilizing effect. It is 
interesting to note that both of the above 
tendencies fluctuate about the final distribution 
at the 210° section. It is to be expected that 


(v’”)§ should be rather strongly affected by curva- 


ture, and measurements of (v’”)3 are necessary for 
further treatment of the problem of curved flow. 


VI. INVESTIGATIONS WITH ACOUSTIC ANALYZER 


The work of Mock and Dryden’ indicates that, 
in general, it would not be possible to express the 
fluctuations as the superposition of a number of 
components of different frequencies. In other 
words, while any single oscillograph record of 
the fluctuations could be represented by a 
Fourier series, it would not be possible to use 
this series to describe another record. An attempt 
to investigate this same problem was made at 


this Institute and gave negative results. An 


acoustic analyzer designed by Mr. L. P. Delsasso 
of the University of California at Los Angeles 
was kindly loaned to the authors. It is essentially 
a quadrant galvanometer with the mirror 
mounted on a trifilar suspension, the funda- 
mental frequency of which could be varied from 
50 to 500 c.p.s. by varying the tension of the 
wires. Since the analyzer is a voltage device it is 
necessary to connect it across a step-up trans- 
former in the output of the amplifier. The hot- 
wire was placed in a fully developed turbulent 
flow and the fluctuations in voltage passed on to 
the amplifier, the output of which actuated the 
analyzer. While a slight response was detected 
at several points in the frequency range attain- 
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Fic. 10. 4/tmean and (u’)*/timean for various positions from fully developed straight to fully 
developed curved flow. For section 9 the abscissa denotes distance from a line parallel to the 
straight portion passing through the center of curvature for the curved portion. For con- 
venience the ordinates for the various profiles are displaced vertically by the amounts indicated. 
Umean = 12.4 meters per second. 


able, the amplitude of the mirror oscillation did 
not reach a steady value, indicating that the 
phase of the disturbances was continually 
changing. Fig. 11 is an oscillograph record of 
fully developed turbulent flow and illustrates the 
apparently random character of the fluctuations. 


VII. Winp TUNNEL MEASUREMENTS 
Experimental determination of the magnitude 
of turbulent fluctuations in wind tunnels is 
highly important because of the fact that certain 
types of aerodynamic measurements are greatly 
influenced by turbulence. For instance Dryden 
and Kuethe’ and later Dryden" showed that in 


%H. L. Dryden, Reduction of Turbulence in Wind 
Tunnels, N. A. C. A. Technical Report 392 (1931). 


fact a marked correlation exists between the 
critical Reynolds number of a sphere and the per- 
centage magnitude of the velocity fluctuations 
present in the air stream. The critical Reynolds 
number was defined as that Reynolds number at 
which the sphere drag coefficient was 0.3. This 
point falls in the vicinity of the center of the more 
or less sudden break which occurs in the curve of 
drag coefficient versus Reynolds number. 

Recent experiments in the GALCIT 10’ tunnel 
have shown that also the maximum value of the 
lift coefficient of airfoils depends largely on the 
degree of turbulence in the air stream. Drs. C. B. 
Millikan and A. L. Klein’ made an extensive 


4C, B. Millikan and A. L. Klein, Effects of Turbulence, 
Aircraft Engineering, August (1933), 
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series of measurements investigating the effect 
of artificially created turbulence on the stalling 
point, or Cronax) Of airfoils. Different degrees of 
turbulence were produced by mounting a grid 
composed of 1/8 inch diameter rods, spaced 3/4 
inch apart, at various distances upstream from 
the model. 

In order to correlate the values of Crmax) 
with the critical Reynolds number of spheres, 
and with the percentage of velocity fluctuations, 
the airfoil model was first replaced by a sphere, 
and the critical drag curves obtained. Then the 
sphere was replaced by the hot-wire, and tur- 
bulence measurements were carried out. Fig. 12 
shows the grid with the sphere in place. Fig. 
13 shows the curves of sphere drag coefficient 
Cp versus R.N. for various distances of the grid 
upstream from the sphere. R.;it, as mentioned 
above, is the R.N. for which Cp=0.3 in each 
case. 

The results of the hot-wire measurements are 
shown in Fig. 14, where (w”)!/u is plotted as a 
function of stream velocity for the clear tunnel 


Fic. 12. The grid in position ahead of a sphere. 
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Fic. 11, Oscillogram of rectilinear fully developed siraight flow taken at the center of the channel. 
wave is shown along with an approximate magnitude scale. tnean = 12.4 meters per second. 
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and two grid positions. In order to correlate 
the values of R-ri: with the turbulence, the values 
of (u’)!/u were taken from the curve of Fig. 14, 
at the same value of tunnel speed at which 
Rerit occurred for the corresponding sphere drag 
test at the same grid position. The correlation 
between R,,;j, and the fluctuations, obtained in 
this way, is shown in Fig. 15. In this figure are 
included the data of Dryden taken at the Bureau 
of Standards. It is remarkable that the values 
fit well on one curve, especially since Dryden 
used several different tunnels, and his apparatus 
was different in detail. 

Fig. 16 shows the correlation between Cymax) 
and the velocity fluctuation (w’”)!/u. The curve 
was obtained in the following way: First curves 
of versus Rerit were plotted for several 
Reynolds numbers of the airfoil. Then the values 
of (u)/u for the different values of Rerit were 
taken from the curve in Fig. 15, and the curves of 


Crunax) versus (u’”)'/u for several R.N.’s of the 
airfoil were plotted. In Fig. 15 the three lower 
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Fic. 13. Sphere drag coefficient drag/(p/2u?-/4D*) 
versus Reynolds number for sphere at various distances 
behind the grid. Reproduced from paper by Millikan and 
Klein, reference 14. 


| 
> 
2, 
het 
os 


A. 
+ 
4 


ia an 


° 20 25 30 40 
vevocity 


Fic. 14. (u”)$/u at various positions behind the grid. 
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Fic. 15. The critical Reynolds number (R-) of a sphere as 
a function of 
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Fic. 16. The maximum lift coefficient, (C1 (max) of an airfoil 
as a function of (u’)3/u. 


black dots indicate only the level of (u”)!/u 
actually measured in the present investigation. 


VIII. FLIGHT MEASUREMENTS 


The discovery of the correlation between the 
maximum lift and turbulence indicated the neces- 
sity of a study of the degree of turbulence en- 
countered by an airfoil in flight. Accordingly, 
Dr. Millikan made an extensive series of free 
air sphere drag measurements in which the 
sphere was mounted on one of the struts of a 
Bird biplane. One might expect that due to tur- 
bulence in the free air, the critical Reynolds 
number of the sphere would be lower than in 
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the wind tunnel. However it was higher than 
that found in the GALCIT wind tunnel and it 
was apparently almost independent of air con- 
ditions and altitude. In view of this fact it 
seemed to be interesting to obtain data on the 
fluctuation of the velocity relative to the plane 
during flight, and to obtain oscillograms of 
atmospheric turbulence. Therefore, a hot-wire 
apparatus was built into a portable form and 
installed in the Bird biplane used for the sphere 
tests. The sphere and hot-wire were mounted 
at corresponding positions with respect to the 
fuselage. 

Due to limitations in space there was little 
opportunity to shield properly the apparatus 
against vibration and electrostatic interference 
from the ignition system, and the hum level 
arising from these factors was so high that it 
was not possible to make accurate velocity 
fluctuation measurements. Oscillograms were 
taken on a special oscillograph designed by 
Mr. L. P. Delsasso and loaned to the authors. 
These oscillograms showed a large percentage of 
periodic fluctuations which almost certainly had 
their origin in either mechanical or electrical 
interference. 

Some estimates were made of the root-mean- 
square value of the fluctuations at various 
speeds, both with power on and in a glide with 
the propeller turning over slowly. Then, for the 
purpose of obtaining a rough approximation to 
the amplifier hum, readings were taken under the 
same conditions with the wire unheated. Sub- 
tracting the hum from the mean square value 
of the output gave a turbulence level approxi- 
mately the same as that in the air stream of the 
GALCIT wind tunnel. However, since it is 
likely that vibrational effects increase with wire 
sensitivity, this value probably represents an 
upper limit of the turbulence. The value found 
is shown in Fig. 15 plotted against the critical 
Reynolds number of the sphere. 


IX. Biimp MEASUREMENTS 


More satisfactory measurements of atmos- 
pheric turbulence were made by installing the 
hot-wire apparatus in the Goodyear blimp Volun- 
teer. The apparatus was well shielded from 
vibrational and electrostatic interference. 
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The mounting for the wire on the blimp con- 
sisted of a 1 inch steel tube 12 feet long lashed 
to the bottom of the control car so that the hot- 
wire, mounted on the end of the tube, was 
situated about 8 feet in front of the control car 
and about 10 feet below the bag. To simplify 
measurement procedure in the air, the hot-wire 
was calibrated immediately before the flight in 
an air stream produced by a variable speed 
vacuum cleaner. 

Unfortunately, the blimp’s stay in Los Angeles 
was very short and it was possible to make only 
one flight, on a day which was quite gusty. 
During the flight it was noticed that the main 
gusts had a period of about 5 to 30 seconds and 
were 3-4 m.p.h. in magnitude. Fluctuations of 
the balancing galvanometer of the Wheatstone 
bridge indicated, in addition, the presence of 
gusts of frequencies between 1 and 10 per second 
and about 1-2 m.p.h. in magnitude. 

The fluctuations to which the amplifier system 
responded (from 50 to 1000 c.p.s.) were very 
irregular. The output meter registered practically 
zero a large part of the time, but during gusts 
made momentary excursions, the maxima of 
which were recorded. The magnitudes of these 
maxima are shown in Fig. 16 where they are 
plotted against the critical Reynolds number of 
the sphere as found in the previous investigation. 
The nature of the higher frequency fluctuations 
is shown in the oscillograph records of which 
Fig. 17 is a sample. For comparison an oscillo- 
gram of the turbulence in the GALCIT wind 
tunnel is shown in the same figure. 

The measurements would seem to indicate 
that velocity fluctuations encountered by air- 
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craft in flight are essentially of low frequency, 
and those of higher frequency, i.e., those which 
are significant for boundary layer stability, are 
of secondary magnitude. 


CONCLUSIONS 


The results obtained, both at this Institute 
and in other laboratories, by means of the hot- 
wire anemometer, indicate the feasibility of this 
method for a study of some of the fundamental 
laws governing turbulent flow. A complete study, 
including flow in various types of channels 
(convergent, divergent, etc.), of (u”)! and, 
especially (v”)!, would add greatly to our knowl- 
edge of this branch of Fluid Mechanics. 

The conclusions drawn, on the basis of the 
results given here, are: 

1. That velocity fluctuations in rectilinear 
fully developed turbulent flow probably obey 


the laws of mechanical similitude as put forth’ 


by von Karman. 

2. That the shearing stress distribution de- 
termines the distribution of the root-mean- 
square value of the velocity fluctuations. 

3. That the fluctuations are of a statistical 
nature and, in fully developed turbulent flow, 
the root-mean-square value is independent of 
channel entrance conditions. 

4. That characteristic curves, such as sphere 


drag or Cramax), can be correlated with (w’”)!/u. 

5. That atmospheric turbulence apparently 
consists of large scale gusts, which themselves 
embody gusts of much smaller scale. The results 
given here, along with the results of Millikan 
and Klein, indicate that only the gusts of smaller 
scale affect boundary layer stability. 


Fic. 17, Oscillograms of wind tunnel and atmospheric turbulence. The upper represents turbulence in the 


free stream of the GALCIT 10’ wind tunnel (V=29 meters per second) (u’*) 1/u=0.4 percent approximate 
time and magnitude scales are shown. The lower record was taken during the blimp investigation (V=19 
meters per second), The tops of a 50 c.p.s, timing wave can be seen at the lower edge. 
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Measurements of the velocity of sound vibrations in ice 
have been made in some detail both by direct and resonance 
methods. The velocity of longitudinal waves was measured 
directly on large bodies of ice. Blasting caps were used to 
create elastic disturbances and the time of travel was 
recorded by an electric seismograph placed at suitable 
distances. A plot of time against distance gives a re- 
markably straight line showing a velocity of 3410 m/sec. 
Laboratory measurements were made on long thin rods of 
ice which were excited by an oscillator and resonance 
frequencies measured both for longitudinal and torsional 


vibrations. The velocities for these vibrations are found to 
be 3163 m/sec. and 1913 m/sec., respectively. A formula is 
derived by which the velocity of longitudinal waves in a 
thin infinite plate can be found from these velocities. This 
formula gives a value of 3400 m/sec. which checks with 
that obtained directly by less than one percent. Various 
other elastic constants for ice are also deduced from these 
measurements. The values of Young’s modulus and 
Poisson’s ratio for ice are found to be 9.17 X10" dynes/cm? 
and 0.365, respectively, in the range —5°C to —15°C. 


INTRODUCTION 


HE measurements made upon the elastic 

constants of ice up to about 1929 gave 
widely discordant results, principally because the 
early measurements were carried out by static 
methods which are ill suited to a substance whose 
elasticity is as imperfect as that of ice. The 
methods which have been applied may be listed 
for convenience in three classes: (a) Static 
methods, including all cases of stretching, bend- 
ing, or twisting in which a steady stress is 
applied and the resulting strain measured; (b) 
Resonance vibration methods, in which the 
natural frequency of a characteristic mode of 
vibration of a small sample of ice is observed; 
(c) Direct velocity determination, in which the 
time for an impulsive elastic disturbance to 
travel a known distance is observed. In the 
present paper there are results obtained by 
the two latter methods. These results agree well 
with each other. 


THE RESONANCE VIBRATION METHOD 


1. Apparatus and method 

The specimens of ice used were formed into 
rods usually a few centimeters in diameter and 
of the order of one meter in length. Both longi- 
tudinal and torsional vibrations were produced 
in these rods. The apparatus used is sketched in 


Fig. 1. For the longitudinal vibrations, a small 
iron disk of mass about 10 g was attached to 
each end of the bar. Short pieces of thin ebonite 
rod frozen into the ends of the bar afforded 
anchorage for the disks. To produce the vibration 
a telephone receiver with the diaphragm re- 
moved was placed near one of the disks and fed 
from the output of a General Radio Low Fre- 
quency Oscillator. A similar telephone receiver 
unit placed near the other disk served to detect 
longitudinal vibrations but was insensitive to 
vibrations of other types. The currents generated 
by this detecting system were fed through an 
amplifier into an a.c. voltmeter. 

For torsional vibrations a small solenoid (mass 
about 15 g) was attached to each end of the 
rod to be tested, being anchored by a small 
piece of ebonite rod frozen into the ice. As is 
shown in Fig. 1 the axis of the solenoid was 


Telephone receiver wut 


from oscillator to amplifier 
and omplifier and volt meter 


Fic. 1. A, side view of apparatus used to obtain longi- 
tudinal vibrations; B, end view of apparatus used to 
obtain torsional vibrations, 
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placed at right angles to the axis of the rod and a 
permanent magnet was placed to produce a 
magnetic field at right angles to both of these 
axes. Thus when alternating current of a suitable 
frequency was fed into the solenoid, torsional 
vibrations of the rod were produced. A similar 
arrangement of solenoid and magnet on the other 
end of the rod was connected through the 
amplifier to the a.c. voltmeter and served to 
detect torsional vibrations while remaining in- 
sensitive to vibrations of other types. 

For both types of vibrations the voltmeter 
which served as detector of resonance gave very 
sharp indications of resonant frequencies. A 
telephone receiver was substituted for voltmeter 
at each resonance setting in order to decide 
whether the rod was resonating to the funda- 
mental note of the oscillator or to an overtone 
of this note. A beat-frequency oscillator was 
found to be very useful for determining the 
approximate frequencies of resonance because of 
the ease with which its frequency could be con- 
tinuously varied. 


2. Results from the resonance vibration method 


Measurements were made upon four rods of 
ice. Rods 1 and 2 were cut from a large block of 
artificial ice, in such a way that the “‘grain”’ of 
the ice—the direction in which the freezing had 
progressed—was perpendicular to the length of 
the rod. Rods 3 and 4 were formed by allowing 
water to freeze slowly and without stirring in 
vertical brass tubes. Since no other precautions 
were taken regarding the crystal orientation, 
the agreement of the results seems to indicate 
that this factor is of little importance in our 
work. The dimensions of the rods were as 
follows: 


Rod Length Cross section 
1 112 cm 5 cm X4 cm 
2 99 cm 6cmX5.5 cm 
3 72.5 cm diam. =3.8 cm 
4 117 cm diam. =3.8 cm 


The results are summarized in Tables I and II. 
The wave-lengths listed are twice the length of 
the rod divided by the number of nodes. The rod 
frequencies, when they differ from the oscillator 
frequencies, are harmonics of the latter. The 
velocities are products of wave-length by fre- 
quency. 
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TABLE I. Velocity of longitudinal vibrations in thin rods of ice, 


Wave- Oscillator Rod Veloc- Resid. 
length frequency frequency ity ual 

Rod cm cycles/sec. cycles/sec. m/sec. m/sec, 
1 224 1386 1386 3105 +58 
1 224 693.1 1886.2 3105 +58 
1 224 467.2 1401 3140 +23 
1 224 354.5 1418.0 3176 —14 
1 224 1387 1387 3107 +56 
1 112 2884 2884 3230 —63 
1 74.7 4231 4231 3161 +2 
1 112 2829 2829 3169 — § 
1 112 1387 2774 3107 +56 
2 198 1593 1593 3154 + 9 
2 198 814.4 1628.8 3225 —62 
2 198 538.7 1616.1 3200 —37 
2 198 407.4 1629.6 3227 —64 
2 99 3255 3255 3223 —60 
2 66 4966 4966 3278 —115 
3 145 2117 2117 3070 +93 
3 145 1100 2200 3190 —27 
3 145 545.5 2182 3164 —- 1 
4 234 1319 1319 3086 +76 
a 234 671 1342 3140 +22 


Ve =3163+9 m/sec. 


TABLE II. Velocity of torsional vibrations in thin rods of ice. 


Wave- Oscillator Rod Veloc- _Resid- 
length frequency frequency ity ual 
cm cycles/sec. cycles/sec. m/sec. m/sec. 
4 234 823.3 823.3 1926 —13 
+ 234 409.2 818.4 1915 —2 
+ 234 275.6 826.8 1935 +21 
4 234 200.9 803.6 1880 +33 
+ 234 163.3 816.5 1911 + 3 


Vr =1914+6 m/sec. 


By use of the well-known formulas 
Vr=(E/p)"?, Vr=(u/p)'” 


we may now obtain all the elastic constants for 
ice from the observed velocities. Thus, taking 
p=0.917 g/cm? 


Young’s modulus, E= 9.17+0.05 x 10" dynes/cm? 
Rigidity, w= 3.3640.02K10"% 
Bulk modulus, k=11.3 +0.6 X10" “ 9 
Poisson’s ratio, o@= 0.365+0.007 

A= 9.1 +0.6 x10" “ 


DirRECT MEASUREMENT OF VELOCITY OF LONGI- 
TUDINAL WAVES IN ICE 


1. Apparatus and method 


The apparatus and method here used were the 
same which have been used extensively in seismic 
prospecting. The waves were generated by an 
electric blasting cap placed in the ice on a lake. 
An electromagnetic seismograph or geophone was 


Bhasting Cap 


Versmograph 


Fic. 2. Schematic diagram of apparatus used in direct 
measurements of velocity in ice. 


placed on the ice in suitable orientation for 
measuring horizontal motion along the line of 
travel of the waves. The currents generated by 
the seismograph were amplified and then re- 
corded in an oscillograph as indicated in Fig. 2. 
The current through the blasting cap was also 
registered on the oscillogram and served to 
determine the instant at which the waves started. 
Timing marks at intervals of 1/100 second were 
placed on the oscillogram by means of a tuning 
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fork. Thus the travel time of the vibrations for 
any desired distance was measured and the 
velocity was determined. The seismograms ob- 


tained at Saylor’s Lake are shown in Fig. 3. 


2. Results of direct measurement of velocity 


One series of measurements was made on 
Saylor’s Lake, a second series on the Lehigh 
Canal. The data obtained are shown in Tables 
III and IV, respectively. The values given in 
these tables for thickness refer to the thickness 
of the ice sheet at the corresponding firing 
station. The readings of time from the oscillo- 
grams were carried out to the nearest thousandth 
of a second but the residual errors are practically 
all less than this amount, indicating a high degree 
of precision for the measurements. The tempera- 
ture of the ice was between 0°C and —5°C in 
each case. Fig. 4 represents the time-distance 
curves obtained from these data. The fact that 
this line is straight and passes through the 
origin is of importance for seismic prospecting. 
In geologic studies by the seismic method the 
time-distance curves obtained are seldom straight 
lines through the origin and it has often been 


TABLE III. Values observed at Saylor's Lake. 


Thickness Distance Time Residual error 


4.5 in. 300 ft. ° 0.027 sec. —0.0003 sec. 
4.5 500 . .044 + 
5.0 700 .061 + .0011 
4.5 900 .080 — .0006 
4.5 1200 .107 — .0004 
4.5 1400 + 
4.2 1600 .142 + 
4.5 1800 .159 + .0009 
4.8 2000 179 — .0013 


Velocity 11,260+20 ft./sec. 


TABLE IV. Values observed at Lehigh Canal. 
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Thickness Distance Time Residual error 
6.0 in 100 ft. 0.009 sec. +0.0002 sec. 
6.0 200 019 — .0005 
6.0 200 .018 + .0005 

300 .028 — 
7.0 400 .038 — .0010 
6.0 500 .046 + .0002 
7.0 600 .055 + .0005 
7.0 700 .065 — ,0003 

700 .063 + .0016 
6.0 800 .075 — .0010 


Velocity 10,800+40 ft./sec. 
Average value of velocity 11,190+30 ft./sec. 
or 3410+9 m/sec. 
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Fic. 3. Seismograms obtained at Saylor’s Lake. 
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Fic. 4. Time distance curve for longitudinal waves in 
ice. Circles, Saylor’s Lake; Crosses, Lehigh Canal. 


argued that some factor such as variation of 
velocity with amplitude of the disturbance was 
responsible. The present work proves that those 
deviations are due to the geologic structure and 
not to some factor inherent in the seismic 
method. 


COMPARISON OF RESULTS 


The velocities Vr, Vp, Vs for longitudinal 
waves in a thin rod, a thin plate, and in an 
infinitely extended solid, respectively, are all 
different. It is well known that 


Our values for Ve and o substituted in this 
formula give 4150+70 m/sec. A value for Vp 
apparently has not been published. The relation 
between the strains €;, €2, €; and the stresses pi, po, 
ps are, for an isotropic solid Ee, = p1—o(p2+ps), 
Ees=ps—a(pitpe). Now, 
for a plane longitudinal wave in a thin plate 
we have p;=0, e2=0, where e; is in the direction 
of propagation of the wave, €2 is in the plane of 
the plate and perpendicular to €), €3 is perpen- 
dicular to the plate. Thus 


or 


Ve=[E/p(1—o*) Ve/ 


Taking the values of Ve=3163+9 m/sec. and 
o = 0.365+0.007 which were obtained from the 
measurements of resonance in thin rods of 
ice the value Vp = 3400+10 m/sec. which checks 
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nicely with the value 
Vp=3410+9 m/sec. 


which is the weighted mean of the values ob- 
tained by direct measurement at Saylor’s Lake 
and Lehigh Canal. 


CORRELATION OF RESULTS WITH OTHER 
OBSERVERS 


Boyle and Sproule! working with longitudinal] 
vibrations of ultrasonic frequencies in thin ice 
bars, obtained a mean value of 3150 m/sec. for 
Vr at 0°C. This value checks exceedingly well 
with our value of 3163 m/sec. at approximately 
—5°C. Their value of 9X10" dynes/cm? for 
Young’s modulus also checks with our value of 
9.17 X10" dynes/cm?. 


In studies on glacier ice by methods similar to 


our direct method, Mothes,? Meisser and Martin, 
and Sorge*® obtained for the velocity of longi- 
tudinal waves an average value of 3500 m/sec. 
This value is very much less than our calculated 
value of longitudinal waves in an infinite solid. 
If glacier ice is less compact than lake ice this 
discrepancy may be accounted for. Also since 
the waves whose velocity was measured, travelled 
quite near the surface, it is possible the ice did 
not behave as an infinite solid. 

K6hler* used a seismic method very similar to 
our direct method for determining the velocity 
of waves in lake ice. He obtained a velocity of 
3.2 k/sec. for longitudinal and 1.7 k/sec. for 
transverse waves in ice 30 cm thick. This value 
for longitudinal waves is in reasonably good 
agreement with that obtained in the present 
work, but Kohler did not take account of the 
fact that the velocity of longitudinal waves is 
different in a thin plate than in an infinite solid. 
This oversight makes his values for the elastic 
moduli of ice unreliable. 
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